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Abstract 

■ We present a discrete total variation calculus in Hamiltonian formalism in 

this paper. Using this discrete variation calculus and generating functions for the 
^ . flows of Hamiltonian systems, we derive two-step symplectic-energy integrators 

Q ' of any finite order for Hamiltonian systems from a variational perspective. The 

, relationship between symplectic integrators derived directly from the Hamil- 

tonian systems and the variationally derived symplectic-energy integrators is 
explored. 

Keywords. Total variation, Hamiltonian formalism, Symplectic-energy in- 
tegrators 
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1 Introduction 

We begin by recalling the ordinary variational principle in Hamiltonian formalism. 
Suppose Q denotes the configuration space with coordinates q l , and T*Q the phase 
space with coordinates (q l ,p l ), i = 1, 2, ■ ■ • , n. Consider a Hamiltonian H : T*Q — > 
| R. The corresponding action functional is denned by 



.s'((./(/).y/(/)Vl : : / (pY-H dt, (1.1) 



Oh 

where [q % (t) , p % (t)) is a C 2 curve in phase space T*Q. 

The variational principle in Hamiltonian formalism seeks the curves (q l (t) , p l (t)) 
r> | for which the action functional S is stationary under variations of (q l (t) , p l (t)) with 

■ fixed endpoints. We first define the variation of (q l (t) , p 1 (t)) . 

Let 

n f> n Pi 



i=l 



be a vector field on T*Q. Here q = (q 1 , ■ ■ ■ , q n ), p — (p 1 , • • • ,p n ). For simplicity, we 
will omit the summation notation YJ in the following. 

Denote the flow of V by F e : F e (q,p) = (q,p), which is written in components as 

q i = f i (e,q,p), (1.3) 
P l =g l (e,q,p), (1.4) 
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where (q,p) £ T*Q and 



d_ 

1c 
d_ 

de 



f{z,q,p) = 4>\q,p)- 



e=0 



e=0 



Let [q l (t) , p l (t)) be a curve in T*Q. The transformations (1.3-1.4) transform 
(q 1 (t) , p l (t j) into a family of curves 

(q\t) 7 f(t)) = (f(e,q(t),p(t)),g*(e,q(t),p(t))). 

Now we are ready to define the variation of (q(t),p(t)): 



8{q\t) lP \t)) 



d 

7k 



(q l (t),f(t)) = (^(q,p),^(q,p)). 



(1.5) 



-=o 



Next, we calculate the variation of S at (q l (t) , p 1 (i)) 
d 



SS 



de 
d_ 
de 
d_ 
~d~c 



S{{q\t),f{t))) 



S(tf\e, q (t),p(t)),g\e, q (t),p(t)))) 



\g\e, q (t), P (t))±f(e,q(t), P (t)) - H(f(e,q(t),p(t)),g\e,q(t),p(t))))dt 



dH 



ij) 1 + -p 



dt+ v 



(1.6) 



If (jf (q(a) , p(a)) = <fi l (q(b) , p(b)) = 0, the requirement of 6S — yields the Hamilton 
equation for (q t (t),p' l (t)) 



^ dp 1 ' 
dH 



(1.7) 



If we drop the requirement of ! (g(a),p(a)) = cb l (q(b) , p(b)) = 0, we can naturally 
obtain the canonical one form on T*Q from the second term in ( |il| . 6 = p^qK 
Furthermore, restricting (q l (t),p' l (t)) to the solution space of ([L7|), we can prove the 
solution of (1.7) preserves the canonical two form u> = dO^ = dp 1 A dq % . 

On the other hand, it is not necessary to restrict (q l (t) , p 1 (t)) to the solution space 
of (1.7). Introducing the Euler-Lagrange one- form 



E(q\f) 

the nilpotency of d leads to 



dp J 



dp* 



dH 



) 



dq 1 



(1.8) 



dE{q\p l ) + -u = 0. 



(1.9) 



Namely, the necessary and sufficient condition for symplectic structure preserving is 
that the Euler-Lagrange one form Q is closed §, g, [kJ]. 

Based on the above variational principle in Hamiltonian formalism and using the 
ideas of discrete Lagrange mechanics (if, |l|, [y| §(], we can develop a natural 
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version of discrete Hamiltonian mechanics with fixed time steps and drive symplectic 
integrators for Hamilton canonical equations from a variational perspective [ flO| . 

However the symplectic integrators obtained in this way are not energy-preserving 
in general because of its fixed time steps [Q, |l8| . An energy-preserving symplectic in- 
tegrator is a more preferable and natural candidate of approximations for conservative 
Hamilton equations since the solution of conservative Hamilton equations is not only 
symplectic but also energy-preserving. To attain this goal, we use variable time steps 
and a discrete total variation calculus developed in [JO] [l2[ [l3[ [l4|, . The basic idea 
is to construct a discrete action functional and then to apply a discrete total variation 
calculus. In this way, we can derive symplectic integrators and their associated energy 
conservation laws. These variationaly derived symplectic integrators are two-step in- 
tegrators. If we take fixed time steps, the resulting integrators are equivalent to the 
symplectic integrators derived directly from the Hamiltonian systems in some special 
cases. 

An outline of this paper is as follows. In Section 2, we present total variation 
for continuous variational principle in Hamiltonian formalism. Section 3 is devoted 
to deriving symplectic-energy integrators. In Section 4, using generating function 
methods, we obtain high order symplectic-energy integrators. We finish this paper 
by making some conclusions and comments in Section 5. 



2 Total variation in Hamiltonian formalism 



In order to discuss total variation in Hamiltonian formalism, we will work with ex- 
tended phase space R x T*Q with coordinates (t, q\p l ) |Q. Here t G R denotes time. 
By total variation, we refer to variations of both (q l ,p l ) and t. Consider a vector field 
on R x T*Q 



d d d 

V = at,q,p)g- t +^(t,q : p)—+^(t, q ,p) — 



(2.1) 



Let F f be the flow of V. For (t,q\p l ) G R x T*Q, we have F e (t,q l ,p l ) = 
which can be written as 



where 



d 
de 
d_ 
de 
d 

7k 



e=0 



t=h(e,t,q,p), 
<f = f (M,g,p), 

P l =5 i (e>*,9,P), 
h{e,t,q,p) = £{t,q,p), 

i 

f(e,t,q,p) = cb l (t,q,p), 

i 

g i {e,t,q,p)=il) i {t,q,p). 



(2.2) 
(2.3) 
(2.4) 

(2.5) 
(2.6) 
(2.7) 



e=0 



The transformation (2.5-2.7) transforms a curve (q l (t) , p l (t)) into a family of curves 
(<f (e,f),f?(e,i)) determined by 



t = h(e,t,q{t),p(t)), 
q i = f(e,t,q(t),p(t)). 
f = g\e.t,q{t),p{t)). 



(2.8) 
(2.9) 
(2.10) 
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Suppose we can solve ( |2.8|) for t: t = h 1 (e,t). Then we have 

q\e,t) =r(e,h-\e 1 t),q(h~ 1 (e 1 t)),p(h- 1 (e,t)))- (2.11) 
=^(e,/i- 1 (e,f),q(^- 1 (e,t)),p(/ l - 1 (6,f))). (2.12) 

Before calculating the variation of S directly, we first consider the first order 
prolongation of V 

prV =£(t, Q,P)^- f + <f>% 9,P)^- + ^'(*, 9>?>)^7 
+ a 4 (t,q,q,p)— +^ 4 (i,q,q,p) — , 

where pr 1 ^ denote the first order prolongation of V and 

a%q,q,p) = D t 4> 1 {t,q,p) - q 1 D t £(t, q,p), 
f3\t, q, q, p) = D t if>% q,p) - fD^t, q, p), 

where D t denotes the total derivative. For example 

D t cf> i (t,q,p) = 4 + <f>iq + (f>i,p. 

For prolongation of vector field and the formula ( 2.14 - 2~15| ), we refer the reader to 



(2.13) 



(2.14) 
(2.15) 



& 



5S = 



Now we calculate the variation of S directly 
d 



de 
d 

d_ 



s((<fM),?M)) 



p*(e,t))^g*( e> t)-fr(g»( e ,t),p»(e,t)) ) dt 



e=0 J a 



e=0 •> a 
6 d 



p l (e, tyj-fk, t) - H@{e, t),p l {e, t)) 



dt 
dt 



■dt 



de 



e=0 
b 



^(e,t))^(e,i)-ff(^(e,t),p'(e,t) ) dt 



(p*(t)q l (t)-H(q l (t),tf(t))) Dt£dt 



b r 



H(q , (t),p t (t)))t+l-p 



dq i J 



dp 1 J 



(i=h(e,t,q(t),p(t))) 



(2.16) 



dt 



+ [p^ l -H{q\p l ))Z]\ a . 
Here in ( 2.16j ), we used ( |2.5| ) and the fact 



(2.17) 



d_ 



e=0 



dt d d 
dt dt de 



t = D t £. 



i=0 



In ( [2.17 ), we used the prolongation formula (2.14). 

If £(a,q(a),p(a)) = Z(b,q(b),p(a)) = and 0*(a, q(a),p(a)) = ^(b, q(b),p(b)) = 
0, the requirement of SS — yields the Hamilton canonical equation 



9" 
P' 



dH 

dp 1 ' 
dH 



(2.18) 
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from the variation (\> l , ip l and the energy conservation law 



(2.19) 



from the variation £. 
Since 



dt 



dH 



dH 

dp 



■P , 



we can easily see that the energy conservation law ( [2.19 ) is a natural consequence of 
the Hamilton canonical equations ( [2.18 ). 
If we drop the requirement 

t(a,q(a),p(a))=t(b,q(b),p(b)) = 0, 
i (o,g(o),p(o))=^(6,g(6),p(6)) = O, 

we ca n define the extended canonical one form on K x T*Q from the second term in 

(HD 



= p i dq l - H(q l ,p l )dt. 



(2.20) 



Furthermor e, res tricting (q l {t) 1 p l {t)) to the solution space of ( 2.18 ), we can prove the 
solution of (2.18|) preserves the extended canonical two form 



dO = dp 1 A dq l - dH{q\p l ) A dt, 
by using the same method in |l5| ]. 



(2.21) 



Remark: The vector field (2.1) can be decomposed along (q l (t) , p 1 (t)) into vertical 
and horizontal components. The vertical component is 

and the horizontal component is 

d d d 

Then analogous to the Lagrangian formalism in Q, the total variation gives rise to 
the Hamilton equations along the vertical direction as well as a relation between the 
Hamilton equations and conservation law for the Hamiltonian along the horizontal 
direction (see also |15[). 
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3 A discrete total variation calculus in Hamiltonian 
formalism and symplectic-energy integrators 

In this section, we develop a discrete version of total variation in Hamiltonian for- 
malism. Using this discrete total variation calculus, we will derive symplectic-energy 
integrators. 
Let 



L{q\p\q\p i )=p i q i -H{q\p i ), 

be a function from K x T(T*Q) to R. Here L does not depend on t explicitly. 

We use P x P for the discrete version of M x T(T*Q). Here P is the discrete 
version of R x T*Q. A point (io)<?0'Po>*i>9i>Pi) ^ P x P corresponds to a tangent 
vector { q }Zf° i ^If" )■ For simplicity, the vector symbols q = (q 1 , • • • , q n ) and p = 

Z\ Zq Z\ Zq 

(p 1 ,--- ,p n ) are used throughout this section. A discrete L is defined to be L : 
PxP^l and the corresponding action to be 



AT-l 

E 

fc=0 



Ht k ,q k ,p k ,tk+i,q k+1 ,p k+1 ){tk+i - t k ). 



(3.1) 



The discrete variational principle in total variation is to extremize S for variations 
of both q k ,p k and t k holding the endpoints (iojQoiPo) an d i^N, Qn'Pn) fixed. This 
discrete variational principle determines a discrete flow <f>: FxP^PxPby 



(3.2) 



Here (ffc+i, Qfc+i,Pfc+i) for all fc G {1, 2, • • • , JV — 1} are found from the following 
discrete Hamilton canonical equation and the discrete energy conservation law 

(tfc+l _ t k)D2Htk,q k >Pk> t k+i,<lk+i>Pk+i) + (*fc - *fe-l)-D5L(tfc_i, 9fc-i:Pfc-i,*fe> 9fc,P fc ) = 0, 
(tfc+i - tfe)-D3L(tfe,q fe ,P fe ,tfe+i,q fe+ i,p fc+ i) + (tfe - tfc_i)£>8L(tfc_i, q fc _i, Pfe_i, tfc, 9 fc ,P fc ) = 0, 

(3.3) 

and 

(tfc+i -*fc)-DiL(tfc,g fe ,p fe ,*fc+i,g fe+ i,P fc+ i) + (** -tfc_i)D4L(tj ! _i,gf fc _ 1 ,p fe _ 1 ,tfc,g fe ,p fe ) 
-L(ife,<7fc,Pfe,ifc+i,<7fe + i,P fc+ i) +L(tk_i,gf fe _ lJ p fc _ 1 ,t fc ,g fc ,p fc ) = 0. 

(3.4) 

Here £)j denotes the partial derivative of L with respect to the ith argument. The 
Eq. ( |3.3| ) is the discrete Hamilton canonical equation that is also called variational 
integrator. The Eq.(3.4) is the discrete energy conservation law associated to([0). 
Unlike the continuous case, the variational integrator (3.3) does not satisfy (3.4) in 
general. Therefore, we need to solve (3.3) and (3.4) simultaneously. 



Now we prove that the discrete flow determined by ( |3.3| ) and (3.4) preserves a 
discrete version of the extended canonical two form to defined in ( 2.21). Therefore, we 
call ( 3.3)-(3.4) a symplectic-energy integrator. We do this directly from the variational 
point of view, consistent with the continuous case p5|. 
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As in continuous case, we calculate dS for variations with varied endpoints. 
dS(t ,q ,p ■ ■ ■ ,t N ,q N ,p N ) ■ (St ,5q ,Sp ■ ■ ■ ,St N ,5q N ,Sp N ) 

N-l 

= ^2 ( D ^ L ( v k) s Qk + D zL(v k )5q k+1 + D 3 L(v k )Sp k + D 6 L(v k )Sp k+1 )(t k+1 - t k ) 



k=0 

N-l N-l 



f J2 (DiL(v k )5t k + D 4 L(v k )6t k+1 )(t k+1 - t k ) + J2 L(v k )(5t k+1 - St k ) 

fc=0 fe=0 
iV-1 

^ (D 2 L(v k )(t k+1 - tk) + D 5 L(v k -i)(t k - t k -i))5q k 

k=l 
N-l 

^2 (D 3 L(v k )(t k+1 - t k ) + D 6 L(v k -i)(tk - tk-i))Sp k 



fe=i 

N-l 

+ ^2 (DiL(vk)(tk+i - tk) + DiL(vk-i)(tk - tk-i) + L(v k -i) - L(v k ))St k 

k=l 

+ D 2 i(«o)(*i - to)Sq + D 3 L(v )(h - t Q )S Po + (DxL(vo)(ti - t ) - L(v ))St 

+ D 5 L(vN-i)(tN - t N -i)Sq N + D 6 L(v N -i)(t N - t N ^i)6p N 

+ (D 4 L(v N ^i)(t N - *jv_i) - L(v N -i))St N , (3.5) 

where v k = (t k ,q k ,p k , t k+ i,q k+1 ,p k+1 ), k = 0, 1, • • • , N — 1. We can see that the 
last six terms in ( |3.5|) come from the boundary variations. Based on the boundary 
variations, we can define two one forms on P x P 

d£(v k ) = D 2 h(v k )(t k+1 - t k )dq k + D 3 h(v k )(t k+1 - t k )dp k 
+ (DiL(« fc )(*fc+i-*fc)-L(«fc))*fc 

and 

9l{v k ) = D 5 h(vk)(t k+1 - t k )dq k+1 + D 6 h(v k )(t k+1 - t k )dp k+1 
+ (D A h(vk){t k +i - t k ) - h(v k ))dtk+i 

Here we have used the notation in We regard the pair (9 7 , Of ) as being the 

discrete version of the extended canonical one form 9 defined in ( |2.20| ) . 

Now we parameterize the solutions of the discrete variational principle by (to, go> ii> Qi), 
and restrict S to that solution space. Then Eq. (3.5) becomes 

dS(t ,q ,p ,--- ,t N ,q N ,p N ) ■ (6t Q ,5q ,Sp - ■ ■ ,8t N ,Sq N ,6p N ) 
= ^(t ,q ,p ,h,q 1 ,p 1 ) ■ (6t a ,Sq Q ,Sp ,St 1 ,Sq 1 ,Sp 1 ) 

+ ^(tN-i,q N -uPN-i^N,qN^PN) ■ ( St N-i,Sq N _ 1 ,Sp N _ 1 ,St N ,Sq N ,Sp N ,) 
= K (^goPo^tfnPi) • (St Q ,Sq ,Sp ,Sti,Sq 1 ,Sp 1 ) 

+ (* Ar_1 )*^L (*o,go:Po>*i>9i)Pi) ' (St Q ,6q ,5p ,5t 1 ,Sq 1 ,Sp 1 ). 

(3.8) 



From (|3.8| ), we can obtain 

dS = 9- + (^ N - 1 r9+. (3.9) 
The Eq. (|1| holds for arbitrary N > 1. Taking N=2 leads to 

<S = 67 + $*6»+. (3.10) 
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Taking exterior differentiation of (|3.10| ) reveals that 
From the definition of 9£ and 9^ , we know that 



- 9+ = dh. 



Taking exterior differentiation of ( 3.12 ), we obtain dO^ — —dO h . Define 

Lu h = dOt = ~d97 . 



(3.11) 



(3.12) 



(3.13) 



Finally, we have shown that the discrete flow <£> preserves the discrete extended canon- 
ical two form wl- 



$*(w L ) = w L . 



(3.14) 



Thus we may call the coupled difference system (3.3)- (3.4) symplectic-energy in- 
tegrator in the sense that it satisfies the discrete energy conservation law ( |3.4| ) and 
preserves the discrete extended canonical two form uj^. 

To illustrate the a bov e discrete total variation calculus, we now present an exam- 
ple. We choose L in (3.1) as 

Htk,q k ,P k ,t k+1 ,q k+1 , Pk+1 ) =p fe+ i^±I Sh. -H(q k+ i, p k+ i), (3.15) 



tk+i — tk 



where p k+ , = 5s±fi±i, q k+ , = 1*±£±L. 



Using (3.3), we can obtain the corresponding discrete Hamilton equation 
Qk+i-Qk-i 1 



P fc+ i -Pk-i 



(tk+i-tk)^(q k+ i,p k+ i) + (tk-t k -i)^{q k _i,p k _i) ) =0, 
(t k+1 -t k ) — {q k+ i,p k+ i) + {t h - t k -i)-^{q k _i,p k _i) ) = 0, 



(3.16) 



where p k _i = Pk+Pk 1 . q k _i = 1 . Using (3.4), we can obtain the correspond- 

ing discrete energy conservation law 

H(q k+ i,p k+ i) = H(q k _ h ,p k t). (3.17) 



The symplectic-energy integrator ( 3.16 )-( 3.17 ) preserves the discrete two form d9^ 
d0 : 



\{dp k A dq k+1 + dp k+1 A dq k ) - H(q k+ i,p k+ i) A 



dt k + dtk+i 



(3.18) 



If we take fixed time steps t k+ i — t k = h (h is a constant), then ( p. 16 ) becomes 



Qk+i - Qk-i 1 ( dH dH \ 

Pk+i - Pk-i 1 (dH dH 



(3.19) 
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Now we explore the relationship between (3.19) and the mid-point integrator for 
the Hamiltonian system 



q = 



p = 



dH 
dp ' 
dH 
dq 



(3.20) 



The mid-point symplectic integrator for (3.2C) is 

<?fc+i -q k _9H 



9fc+i q k 



dp 
dH 
dq 



(q k+ i,p k+ i) 



(3.21) 



In ( |3.21 ), we replace k by k — 1 and obtain 

q k - q fc _i dH 



Pk-Pk-i 



^-(q fe _i,p fc _i), 
dH 

-^-(9fe-i,Pfc-i)- 



(3.22) 



Adding ( 3.22| ) to (3.21) results in (3.19). Therefore, if we use ( [3.21 ) to obtain q k , 
the two-step integrator (3.19) is equivalent to the mid-point integrator (3.21). How- 
ever, the equivalence does not hold in general. For example, choose L in (3.1) as 



Htk,q k ,p k ,t k+1 ,q k+1 ,p k+1 ) = p k 



<?fc+i - qk 



tk-\ 



tk 



H(q k+ i,p k+ i] 



(3.23) 



and take fixed time steps t k +i —t k ~h. Then (3.3) becomes 

Qk+i - qk i 



dH 

dp 

dH 



Pk - Pfc-i 



dq 



(q k+ i,p k+ i 



dH 
dp 
dH 



dq 



(q k -l,Pk-l) 



(3.24) 



The integrator ( 3.24 ) is a two-step integrator which preserves the two form dp k Adq k+ i . 
In this case, we cannot find one-step integrator that is equivalent to ( 3.24 ). 

In conclusion, using discrete total variation calculus, we derive two-step symplectic- 
energy integrators. When taking fixed time steps, some of them are equivalent to 
one-step integrator derived directly from the Hamiltonian system while the others do 
not have this equivalence. 



4 High order symplectic-energy integrators by gen- 
erating functions 

In this section, we develop high order symplectic-energy integrators by using the 
generating function of the flow of the Hamiltonian system 

z = J V H(z), (4.1) 



where z = (p, q) T , J = 
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Wc first recall the generating function with normal Darboux matrix of a symplectic 
transformation . For details, see || 0]. 

Suppose a is a An x 4n nonsingular matrix with the form 

A B 
C D 

where A, B, C and D are both 2n x 2n matrices. 
We denote the inverse of a by 



a- 1 



Ai Bi 
Ci Dt 



where Ai,B±, C\ and D\ are both 2n x 2n matrices. 
We call a 4n x 4n matrix a a Darboux matrix if 

a T J4„o; = J4„, (4.2) 

where 

, -I 2n \ f (Jin \ T f -I T , 
■Jin — r n I 5 '-'4? 



/an y ' -™ v - J 2«; ' z " \In 

where /„ is an n x n identity matrix and /2n is a 2n x 2n identity matrix. 

Every Darboux matrix induces a fractional transform between symplectic and 
symmetric matrices 

<7 Q : Sp(2n) — > Sm(2n), 

<r a (S) = (AS + B)(CS + D)- 1 = M, for 5 G Sp(2n), det(CS + D) ^ 

with the inverse transform cr^ 1 = cr ct _ 1 

o-^ 1 : 5m(2n) -> Sp(2n), 

a- x (M) = (AiM + Bi)(CiM + Dx)" 1 = 5, 

where Sp(2n) is the group of symplectic matrices and Sm(2n) the set of symmetric 
matrices. 

We can generalize the above discussions to generally nonlinear transformations on 
K 2li . Denote Spnl(2n) the set of symplectic transformations on M. 2n and Smnl(2n) 
the set of symmetric transformations, i.e. transformations with symmetric Jacobian, 
on R 2n . Every / £ Smnl(2n) corresponds, at least locally, to a real function </> (unique 
to a constant) such that / is the gradient of <p: 

f(w) = V 0(u>), (4.3) 

where y<t>(w) = (<f> Wl (to), • • ■ , </>™ 2 „ (u>))< 
Then we have 

<7 Q : Spnl(2n) — > Smnl(2n), 

a a (g) — (A o g + B) o (C o g + D)^ 1 = V 0, 

for g £ Spnl(2n), det(C<? z + Z?) ^ 0, or alternatively 

A 5 (*) + J B*=( V 0)(C 5 (*)+23a!) 1 
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where o denotes the composition of transformation and the 2n x 2n constant matrices 
A, B, C and D are regarded as linear transformations. g z denotes the Jacobian of the 
symplectic transformation g. 

We call 4> the generating function of Darboux type a for the symplectic transfor- 
mation g. 

Conversely, we have 

ct" 1 : Smnl(2n) — > Spnl(2n), 

o"a X (V0) = (4l ° + Si) o (d o V + D^- 1 = g, 

for det(Ci</) u , l „ + D x ) ^ 0, or alternatively 

A x V + B x w = g{Ci v + D x w), 

where g is called the symplectic transformation of Darboux type a for the generating 
function cj). 

For the study of integrators, we may restrict ourselves to the normal Darboux 
matrices, i.e., those satisfying A + B = 0, C + D = I% n , The normal Darboux matrices 
can be characterized as 

J E h n J -E)> E =\^ + J ^F), F T = F, (4.4) 
and 

-1 _ ( (E - hn)J2n hn\ ^ ^ 



(X 



The fractional transform induced by a normal Darboux matrix establishes a one- 
one correspondence between symplectic transformations near identity and symmetric 
transformations near nullity. 

For simplicity, we take F = 0, then E = hl% n and 



J2n —J; 

1 T 1 T 



(4.6) 



Now we consider the generating function of the flow of (4.1). Denote the flow of 
(4.1) by e l H . The generating function <p(w,t) for the flow e l H of Darboux type (4.6) 
is give by 

= {'hn »4 - J2n) ° (\ e \l + \ l2 ^j , for small \t\. (4.7) 
<p(w,t) satisfies the Hamilton- Jacobi equation 

^{w,t) = -H (w + ^J 2n ycj>{w,t)\ (4.8) 
and can be expressed by Taylor series in t 

oo 

4>{w,t) = ^cj) k (w)t k , for small \t\. (4.9) 
fc=i 
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The coefficients (f> k (w) can be determined recursively 
<f>\w) = -H(w), 



m=l jiH hj m =fc 



where /c 1, and we use the notation of the m-linear form 



From (4.7), we can see that the phase flow £ := e l H z satisfies 



oo 

J 2n (i-^) = V^( £ T £ ) = E* i V^'(^ l£ )- (4.11) 



Now we choose L in (3.1) as 



L(t fc ,q k , Pk ,t k+1 , q k+1 , p k+1 )= p k+ x ^ - (<Z fe+ 1 , P fe+ 1 ) , (4.12) 



where 



V/"(g fe+ i, Pfe+ i) = XV^(9*+*>P*+i)- (4.13) 

3=1 



The corresponding symplectic-energy integrator is 
2 2 ( v ( * fe + 1 _ + (** - **-0-^-(9*-i,P*-i) J = 0, 

V> m (<7 fe+ i ,p fc+ i) = ^ ra (q fc -i,p fc -i), 

(4.14) 

which preserves the discrete extended canonical two form 

^dp k Adq k+1 +dp k+1 Adq k )-r\q k+ i,P k+ ^^(^^^) ■ (4.15) 



The integrator ( 4.14 ) is a two-step symplectic-energy integrator with 2mth order of 
accuracy. 

5 Concluding remarks 

We have developed a discrete total variation calculus in Hamiltonian formalism in this 
paper. This calculus provides a new method for constructing structure-preserving 
integrators for Hamiltonian systems from a variational view of point. Using this 
calculus, we have derived discrete energy conservation laws associated to the integra- 
tors with variable time-steps. The coupled integrators are two-step integrators and 
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preserve a discrete version of the extended canonical two form. If we take fixed time- 
steps, the resulting integrators are equivalent to the symplectic integrators derived 
directly from the Hamiltonian systems only in special cases. Thus, some new two-step 
symplectic integrators have variationally been obtained. Using generating function 
method, we have also obtained high order symplectic-energy integrators. 

In principle, our discussions can be generalized to multisymplectic Hamiltonian 
system 



Mz t + K z x = VzH(z), zeR n , 



(5.1) 



where M and K are skew-symmetric matrices on R n ,n > 3 and S : R™ — > R is 
a smooth function |^|. We call the above system multi-symplectic Hamiltonian 
system, since it has a multi-symplectic conservation law 



d 



d 



dt dx 

where lo and k are the pre-symplectic forms 

lo = —dz A Mdz, 
2 

Constructing action functional 

'1 



K = 0, 



k = —dz A Kdz. 
2 



(5.2) 



S = 



-z 1 (Mz t + Kz x ) - H(z) dx Adt 



(5.3) 



and performing total variation on (5.3), we obtain the multisymplectic Hamiltonian 
system (5.1) as well as the corresponding local energy conservation law 



|(s ( .)-i^) + A(i^,_o, 

and the corresponding local momentum conservation law 



d_ 

dt 



l 



-z 1 Mz x + 



d_ 

d.r 



S{z) - X -z T Mz t 



0. 



(5.4) 



(5.5) 



In the same way, we can develop a discrete total variation in multisymplectic for- 
mulation and obtain multisymplectic-energy-momentum integrators. We shall explore 
this issue elsewhere. 

On the other hand, we may also get symplectic/multisymplectic-energy-momentum 
integrators by means of the difference variational principle || , |U)[ in view of regard- 
ing difference with variable steps as an entire object. In this approach, the difference 
Legendre transformation may be introduced so as to the discrete total variations in 
Lagrangian and Hamiltonian formalism may be transformed to each other. We shall 
also explore this issue elsewhere. 
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